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Applying chiral perturbation to twisted mass lattice QCD*^ 

Stephen R. Sharpe a and Jackson M. S. Wu a 

a Physics Department, Box 351560, University of Washington, Seattle, Washington 98195, USA 

We have explored twisted mass LQCD (tmLQCD) analytically using chiral perturbation theory, including 
discretization effects up to 0(a ), and working at next-to-leading (NLO) order in the chiral expansion. In 
particular we have studied the vacuum structure, and calculated the dependence of pion masses and decay 
constants on the quark mass, twisting angle and lattice spacing. We give explicit examples for quantities that 
both are and are not automatically improved at maximal twisting. 



1. Introduction 

There has been much interest recently in the 
twisted mass formulation of lattice QCD I1I2| (see 
Ref. 3 for a recent review). We study here the 
importance of the symmetry breaking inherent in 
tmLQCD, which can be done analytically using 
the chiral effective theory including effects of dis- 
cretization. 

2. The effective chiral Lagrangian 

To derive an effective continuum chiral theory 
for tmLQCD, we follow the two-step procedure of 
Ref. 4 . We first write down an effective contin- 
uum Lagrangian at the quark level that describes 
the long distance physics of the underlying lattice 
theory, and then match it onto an effective chi- 
ral Lagrangian. The details of this procedure are 
given in Ref. p]. 

The effective continuum Lagrangian is con- 
strained to be invariant under the symmetries of 
the lattice theory. It has form 3 : 



+ bxa^ia^F^ifj + 0(a 2 ) , 



(1) 



where C g is the continuum gluon Lagrangian, m 
is the physical quark mass, defined as usual by 

'Partially supported by the U.S. Department of Energy 
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tTalk presented by J. Wu at Lattice 2004, Fermi National 
Accelerator Laboratory, June 21 - 26, 2004. 
3 The 0(a 2 ) terms do not contribute to any further sym- 
metry breaking than the terms explicity shown. 



m = Z m (mo — fh c )/a, and [i is the physical 
twisted mass defined by \i = Z^/io/a (mo and 
Ho are bare normal and twisted mass parameters 
respectively). Note that the lattice symmetries 
forbid additive renormalization of /iq p^. 

In writing down JQ), we have dropped terms 
that vanish by equations of motion. We have also 
dropped terms of order higher than quadratic in 
{to, [i, a} (factors of Aqcd are implicit here and 
henceforth) in anticipation of the power-counting 
in our chiral effective theory. The net result is 
that 0]) differs from that for the standard Wil- 
son theory given in 0] only by the addition of a 
twisted mass term. 

Next we write down a generalization of the con- 
tinuum chiral Lagrangian that includes the effects 
of the Pauli term. We use the power counting 
scheme: 

1 > {m,p 2 ,a} > {m 2 , mp 2 ,p 4 , am, ap 2 , a 2 } , (2) 

where m is a generic mass parameter that can be 
cither m or [i. The chiral Lagrangian is built from 
the standard SU (2) matrix- valued field E. It can 
be obtained from the quark Lagrangian, Q), by 
a standard spurion analysis. We must introduce 
a spurion matrix A for the Pauli term, as well 
as the usual spurion x for the mass terms. Both 
transform in the same way as the S field. At the 
end of the analysis the spurions are set to their 
respective constant values: 



X 

A -» 2ir a 



2£>q(to + ir 3 /i) = m + ir 3 jj, , 



(3) 



1 



2 



where B and Wo are unknown dimensionful pa- 
rameters. The difference from the standard con- 
struction comes only from the inclusion of the 
/i term in the constant value of x- Thus, we 
can read off the form of the chiral Lagrangian 
for tmLQCD from Ref. |Hj , where the Lagrangian 
for untwisted Wilson fermions was worked out to 
quadratic order in our expansion 4 . The only ex- 
tension we make is to include sources for currents 
and densities. 

3. The phase diagram of tmLQCD 

Using the chiral Lagrangian, we have extended 
the analysis of Ref. 0] into the twisted mass 
plane. The details of the analysis are given in 
Ref. U3- We record here the main results. 5 

In the region a 3> to' ~ /i ~ a 2 where the 
Aoki phase occurs for untwisted Wilson fermions 
(to' = to + a Wq/Bq is the shifted mass due to 
effect of the Pauli term, as noted in Ref. 0]), the 
potential energy is dominated by: 



C3/C2 ~ n/a 2 



V x = -^Tr(£ + £t) 

, c 2 



|Tr[i(E-St) T3 ] 



^[Tr( S + £t )r , 



(4) 



where c\ ~ to', C3 ~ /i, and C2 ~ a 2 are all of 
the same order. As in the untwisted case, it is 
the competition between the leading order mass 
terms and the a 2 term at NLO that can lead to 
interesting phase structure. 

To determine the vacuum structure we find the 
condensate, £01 which minimizes the potential. 
Parametrizing using T.q — A m + iB m • r, where 
A 2 ^ + B^ n = 1, the functional forms of A m and 
B m are found by solving a quartic equation. The 
minimum is found to lie in the I-T3 plane in flavor 
space. 

The vacuum phase structure depends on the 
sign of the coefficient C2, which is a linear com- 
bination of low-energy constants (LEC's) in the 
chiral Lagrangian. Fig. ^ illustrates the two pos- 
sibilities. 

For C2 > 0, the usual Aoki phase (first noted in 
Ref. U|) occurs, but is confined to a short segment 

4 Similar work has been done in Ref. |7|. We have extended 

that work by including 0(a 2 ) terms. 

5 Similar considerations have been made in Refs. 1811 II 
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(a) Phase diagram for C2 < 
C3/C2 ~ n/a 2 
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(b) Phase diagram for C2 < 



Figure 1. Phase diagram of tmLQCD. The sign 
of the coefficient c-i determines whether flavor 
symmetry is spontaneously broken in the stan- 
dard Wilson theory. The solid lines are first order 
transitions across which the condensate is discon- 
tinuous, with second-order endpoints. 



on the untwisted axis. In the Aoki phase segment, 
the charged pions are massless but the neutral 
pion is not. All three pions become massless only 
at the second-order endpoints. 

For C2 < 0, there is a first-order transition line 
extending a distance A/j, ~ a 2 in the twisted di- 
rection. The charged pions are not massless any- 
where on this first-order line (including the end- 
points), while the neutral pion becomes massless 
at the two second-order endpoints. In fact the 
neutral pion mass for ci < has the same depe- 
nence on \i along the first-order transition line as 
that on to along the Aoki phase line for ci > 0. 
Evidence for the phase structure in the ci < 
scenario has recently been seen in numerical sim- 
ulations |1U) . 



3 



4. Pion masses and pion decay constant 

In the region m^fi^a^a 2 ^ am ~ m 2 
away from the "Aoki phase" (where most simu- 
lations will be done), we find that the splitting 
of pion masses due to parity-flavor breaking first 
occurs at NLO: 



a 2 2 2 C d j.1 

Am, = m„ — = ^ 

Wl - 3 (m + a) 2 +ft 2 



(5) 



where C is a linear combination of the LEC's. The 
mass splitting is oc a 2 for any values of m and ft 
within the region we are considering, not just at 
maximal twisting. For ffi,/i > a, (JSJl becomes 
Ca 2 sm 2 u>, the form expected from the symme- 
tries of tmLQCD 0. Note that at NLO, there are 
only tree level contributions to the mass splitting, 
whose full expression agrees with that of Ref . . 
For pion decay constants, we find no splitting un- 
til next-to-next-to-leading order. 

The results of the one-loop calculation of the 
pions masses, the wavefunction renormalizations, 
and the pion decay constants will be given else- 
where |12| . 

5. Current matrix elements 

Once the twisting angle is defined, the physical 
untwisted currents can be obtained simply from 
the appropriate functional derivative of the effec- 
tive chiral Lagrangian with respect to the sources. 
We use the non-perturbative definition given by: 



tanc = (dtfP 1 )/(d ll A] l P 1 ), 



(0) 



where the defining correlators are evaluated in the 
chiral theory using single pion contributions. We 
find the twisting angle is related to the parame- 
ters in the chiral theory by: 



tanw = B 3 , m /(A m + 4W w a/f 2 ) 



(7) 



where Wio is one of the unknown LEC's. 

As an example, consider the physical untwisted 
axial vector current: 



[A$(l>) = (co S (u;)(l-5 3a ) + 5 3a )A« 



+ sin(w)e 3Qb y b . (8) 



The effect of parity-flavor breaking can be seen 
when [A a ^\ is expanded in terms of the pion fields. 
There are terms with even number of pion fields 
that are usually forbidden by parity, and the co- 
efficients of the expansion depend on the flavor 
index. 

The physical, single pion matrix element of 
the axial vector current, (0| [j4"] |7Tb), which de- 
fines the pion decay constant, is an example of 
a quantity that is automatically 0{a) improved. 
But the unphysical, two-pion matrix elements are 
not (and are not expected to be, according to 
the general theory of Ref. For instance, let 
q = p2 — pi <C m 2 ,, then for a = 1, 2, wc find: 



(n a (pi)\[A a u ]\n 3 (p2)) = -Opi 



■P2) — 



+iq 



4B 



v 2 



SB. 



3,m 



r- 

A3 



aV 3 



tVa 



(9) 



p p 

where the T>i's are linear combinations of LEC's. 
The 1 /m 2 . term comes from a diagram which con- 
tains a three-pion vertex that arises from parity- 
flavor breaking. Since the continuum limit must 
be taken before the chiral limit to avoid entering 
the Aoki phase, there is no divergence. 

A potential use of the unphysical two-pion cur- 
rent matrix elements is to provide a way of calcu- 
lating the new LEC's introduced by discretization 
into the chiral Lagrangian. 
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